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We study the T(35, 25") -* rjT{lS) and T(3S, 25") -* tt + tv-Y(1S) transitions with 122x 10 6 T(35) 
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PACS numbers: 14.40.Pq,13.25.Gv 

The QCD multipole expansion (QCDME) model [J 
describes hadronic transitions between heavy quarkonia. 
Despite its success for hadronic transitions in charmo- 
nium, this model has limits in explaining all hadronic 
transitions in the bottomonium spectrum. The QCDME 
predicts the suppression of the transitions between bot- 
tomonia via a rj meson with respect to those via a dipion, 
the former being associated with the spin-flip effects of 
the b quark. The Y(4S) -> rflT(lS) and Y(2S) -> rflT(lS) 
transitions have been observed at rates significantly dif- 
ferent from the predicted values 0, The measured 
width T[Y(2S) — > r]T(lS)] is smaller than predicted, 
while T[Y(4S) -> rjT(lS)} is larger than T[T{4S) 
7r + 7r _ T(15)], although it was expected to be suppressed 
in analogy with decays of the lower-mass T resonances. 
Some suggest that the latter result could be related to 



above-SB threshold effects 0,0. The T(3S) -> rjY(lS) 
transitions have not been observed 0. Precise measure- 
ments of the transitions between bottomonia via a r] me- 
son, as well as their rate with respect to the dipion tran- 
sitions could shed light on the chromomagnetic moment 
of the b quark. 

In this paper we study the transitions Y(nS) — > 
T]Y(1S) and Y(nS) -> 7r+7r-T(15) with n = 3, 2 
and measure the ratios of partial widths T[Y(nS) -^>- 
rfT(lS)]/T[T(nS) -> w+tt-Y(1S)]. The transitions are 
studied for events in which the Y(1S) decays to either 
or e + e~. The r\ meson is reconstructed from its 
77 and 7r + 7r~7r° decay modes, where the 7r° decays to 77. 
The analysis thus considers the final states 7T 4 tt~ '77^ t~ , 
"f"f£ + £~ and 7t + tt~£ + £~ , where £ = e or fi. 

We analyze BABAR data samples consisting of (121.8 ± 
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1.2) x 10 6 Y(3S) and (98.6 ± 0.9) x 10 6 Y(2S) mesons. 
These correspond to integrated luminosities of 28.0 fb _1 
and 13.6 fb _1 , respectively. We use 2.6 fb _1 collected 
30 MeV below the Y(3S) resonance, and 1.4 fb _1 col- 
lected 30 MeV below the Y(2S) resonance ("off-peak" 
samples) for background studies. 

The BABAR detector is described in detail elsewhere [H, 
. We briefly mention the features relevant to this anal- 
ysis. Charged-particle momenta are measured in a five- 
layer double-sided silicon vertex tracker (SVT) and a 40- 
layer central drift chamber (DCH), both embedded in 
a 1.5-T axial magnetic field. Charged-particle identifi- 
cation is based on specific energy loss in the SVT and 
the DCH and on measurements of the photons produced 
in the fused-silica bars of a ring-imaging Cherenkov de- 
tector. A CsI(Tl) electromagnetic calorimeter (EMC) is 
used to detect and identify photons and to identify elec- 
trons, while muons are identified in the instrumented flux 
return of the magnet (IFR) . 

Monte Carlo (MC) simulated events, used for efficiency 
determination and selection optimization, are generated 
with EvtGen [§]; GEANT 4 [§] is used to simulate the de- 
tector response. The variations of conditions and beam 
backgrounds are taken into account in the simulation. 
The simulated events are then analyzed in the same man- 
ner as data. Large MC samples simulating inclusive 
Y(3S) and T(2S) decays including all known and pre- 
dicted transitions and continuum e + e~ — > e + e~(-/) and 
e + e~ — > /i + /i _ (7) processes are used to characterize the 
backgrounds. Background from continuum quark pro- 
duction is negligible. In the MC signal samples, the dis- 
tribution of generated dilepton decays incorporates the 
Y(IS) polarization. Dipion transitions are modeled ac- 
cording to the matrix elements measured by CLEO ficj . 
The angular distribution in Y(3S, 2S) — > riY(lS) pro- 
cesses is generated as a vector decaying to a pseudoscalar 
and a vector. The r\ — > 7r + 7r~7r° decays are modeled ac- 
cording to the known Dalitz plot parameters [ll|. Final 
state radiation effects are described by PHOTOS [lgj . 

Events of interest contain two oppositely charged par- 
ticles, identified as either electrons or muons. A fit con- 
strains them to originate from a common vertex and 
to have invariant mass Mgg equal to the known Y(1S) 
mass [llj. The fit must yield a x 2 probability > 10~ 5 . 
Muon identification is based on the energy deposited 
in the EMC, and the number of coordinates and in- 
teraction lengths traversed in the IFR. Electron iden- 
tification is based on specific energy loss in the SVT 
and DCH combined with energy deposition in the EMC. 
Bremsstrahlung energy loss is partially recovered by an 
algorithm combining the energy of an electron candidate 
with the energies of nearby photons. 

Besides the lepton pair, we require a pair of oppositely 
charged tracks not identified as electrons and/or two neu- 
tral particles identified as photon candidates. Events 
with additional charged tracks are rejected. A fit con- 



strains all final state particles to originate from a com- 
mon vertex, to have a total energy equal to the sum of the 
beam energies, and an invariant mass equal to the Y(3S) 
or Y(2S) mass [Uj. The fit must yield a \ 2 probability 
> 10~ 5 . 

A trigger-level prescaling of Bhabha scattering events, 
whose signature is given by two electrons of large invari- 
ant mass and no additional charged track of transverse 
momentum > 250MeV/c, causes the efficiency for the 
final states containing electrons to be smaller than for 
final states with muons. The di-electron efficiency drops 
to ~ for Y(2S) transitions in all final states considered, 
and it is < 3% for Y{3S) -)• 7)Y(1S) in the 776+^ final 
state. These final states are not considered further. 

The event selection criteria have been optimized sep- 
arately for each final state. The background contribu- 
tions have been studied using MC samples of inclusive 
Y(nS) decays and of e + e~(7) and ^+^^(7) events. The 
MC background yield has been compared to real back- 
ground yield from data and found to be compatible with 
it within the uncertainties. Also it has been verified that 
the distributions of all the discriminating variables are 
well-described by the MC background. 

No further selection is applied for the ir + ir~£ + £~ final 
states, while the additional requirements summarized in 
TableUare needed for the other final states. To select the 
7r + 7r~77i' + i'~ final states we require that the two-photon 
invariant mass M 77 be compatible with the ir° mass. 
Background events are rejected by applying selection cri- 
teria to the opening angle between the two pions, calcu- 
lated in the e + e~ center-of-mass (CM) frame (0*^), and 
also to the invariant mass of the dipion candidate calcu- 
lated assuming the electron mass hypothesis (M conv ). In 
particular for the 7r + 7r _ 77e + e _ states, M conv > 30 MeV 
suppresses events in which a photon converts in the de- 
tector material and the electrons are reconstructed as pi- 
ons. In the T(3S') — > t]Y(lS) final state, cross-feed from 
Y(3S,2S) — > iT + Tr~Y(2S,lS) transitions is suppressed 
by vetoing events with AM^n = M^nii — Mu compatible 
with any of the known mass differences between narrow Y 
resonances. In the 77/i + /i~ final states, the backgrounds 
due to the radiative transitions Y{nS) — > 7X(,j(2P, IP) 
with Xftj(2P, IP) — > jY^S) are rejected by vetoing 
events where either photon energy calculated in the CM 
frame (E* fl 72 ) is compatible with any of those transi- 
tions. The background from /i + /i~(7) events is reduced 
by requirements on the opening angle between the two 
photons (6** 7 ), and on the momentum of each lepton 
(p| ± ) in the CM frame. 

The signal yields are extracted with a 2-dimcnsional, 
unbinned, extended, maximum-likelihood fit to the mea- 
sured distribution of a pair of variables. For the Y(nS) — > 
tt + tt~Y(1S) transitions we fit the AM„ versus Mu 
distribution, both calculated from the measured values 
prior to the invariant mass constraint (Fig. [TJ). For the 
Y(nS) —> T]Y(1S) transitions with 77 — > ir + ir~ir° de- 
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TABLE I: Additional 

■K + n~"f"ft r l" and yy£ + £ 



requirements applied to select 
final states. Masses are expressed 
in MeV/c 2 , energies in GeV, and momenta in GeV/c. 



T(2S) -> yTjlS) 



7T 1 TV 77^'^ 



7T + 7T 77/i ' /i 



r(3g) 7?r(ig) 



7r 7r 77e e 



M, 



< 310 



Mr 



< 280 



30 < M conv < 280 



90 < M 77 < 180 



cos 0* < 



100 < M 77 < 170 90 < M 77 < 150 



(400 < AM„ < 550) U (AM„ > 580) 



77/7^ 



200 < £ 7l , 72 < 350 (120 < E*]_ < 360) U (490 < E* fl < 660) 
(130 < E* 2 < 370) U (470 < E* 2 < 700) 

cos 8* < —0.86 



4.60 < p* t± < 4.85 



-0.45 < cos6> 77 < 0.22 
4.4 < p* + < 5.0; 4.4 < p*_ < 5.1 



cays, we fit the AM n versus M W7r77 distribution, where 
AM, = M Tir77 « — Mm — Mnn-yy and M 7r7r77 is the in- 
variant mass of the r/ decay products (Fig. [2]). For 
the T(nS) — > r)T(lS) transitions with -q — > 77 de- 
cays, we fit the AM^ versus M 77 distribution, where 
AM^ = M 77 « - M/, 



M 77 (Fig. El). 



Each observed distribution is fit to a sum of a signal 
and a background component. The functional form of the 
probability density functions (PDFs) for signal and back- 
ground have been determined from MC samples. The 
signal PDFs are described by double or triple Gaussian 
functions, or by a Gaussian-like analytical function with 
mean value /i but different widths, <jl,r, on the left side 
(for x < n) and on the right side (for x > fi) plus asym- 
metric tails oll,r, defined as: 



J-(x) = cxp I — 



(x /i) 2 



2a L,R + a L,R(x - A 1 ) 2 



(1) 



The PDFs used to model the signal and background 
shapes in each fit are given in Table ITT1 The free parame- 
ters in each fit are the signal and background yields and 
the parameter of the background PDFs of Table [TTJ The 
signal shape parameters are also floated in the fits to the 
T(nS) —> tt + tt~£ + £~ samples, while they are fixed to the 
values determined from MC samples in all other cases. 

The number of signal candidates returned by the fits 
is reported in Table IIII1 We estimate the signal sig- 
nificance in standard deviations as */ 2 log[£(AT) / C(0)} , 
where £(N) / C(0) is the ratio between the likelihood val- 
ues for a fit that includes a signal yield N and a fit 
with a background hypothesis only. For the T(2S) — > 
rpT(lS) transition the signal significance is 8.0<r for the 
7r + 7r~77/! + /!~ final states and 8.5<r for the 77/i + /i _ ones. 
For the Y(3S) — > ijY(lS) we find no evidence of a signal 
in any of the final states considered and calculate 90% 
CL upper limits on the number of signal events (Njjl) 
as f* UL C(N)dN = 0.9 x J™ C(N)dN. The efficien- 



cies with which signal events satisfy the selection crite- 
ria (e se i) are determined using simulated signal samples. 
Corrections are applied to account for differences between 
data and MC in lepton identification and 7r° reconstruc- 
tion efficiencies. The corrected values are also reported 
in Table [TTTJ 



s, 



5000- 




0.56 



0.58 



AM^CCeY/./) 6 





0.9 0.92 
AM^GeV/c 2 ) 




0.9 0.92 
AM^GeV/c 2 ) 




9.4 9.6 
M ee (GeV/c 2 ) 



FIG. 1: AM„ and M u distributions for (a,b) Y(2S) -s> 
■k + -k~T(1S) — > ir + tt~ n + /i~ candidates, (c,d) T(3S) — > 
tv + tt~T(1S) — s> 7r + 7r~/i + /i~ candidates, and (e,f) T(3S) — > 
■k + -k~T(1S) — > -K + -K~e + e~ candidates. Data are represented 
by dots, the fit results as solid curves and the background 
components by the dashed curves. 

Possible sources of systematic uncertainty are con- 
sidered in addition to those on the number of Y(nS) 
(N?) and on the values for secondary branching fractions 
(B seC ondary) • The uncertainties on charged-particle 
track and single 7 or tt° reconstruction efficiencies are de- 
termined by a comparison between data and MC events 
using independent control samples of r pair events, each 
r decaying to either one or three charged-particle tracks. 
The systematic uncertainty on the muon or electron iden- 
tification probability is estimated by comparing the val- 
ues determined in the Y(nS) — > 7r + 7r _ T(15 l ) mode in 
data and MC samples. For each discriminating variable, 
we compare the distribution for the signal component 
deconvolved from data with the maximum likelihood fit 
used for the extraction of the yields [3] to the distri- 
bution obtained in the MC. The related systematic un- 
certainty is estimated as the change in event selection 
efficiency induced by the difference between the distribu- 
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TABLE II: Functions used to model the signal and background PDFs. 



Final state 


1 st Variable 


Signal 


Background 


2 nd Variable 


Signal 


Background 


T{nS) - 


->■ Tv + ir-£ + r 


AM„ 


triple Gaussian 


th order poly 


M« 


Eq.© 


th order poly 


r(2S) - 


■¥ 7r + 7r _ 77/i + /i _ 


AM, 


triple Gaussian 


th order poly 




Eq.© 


l st order poly 


r(2S) - 


•> 77M + M~ 


am; 


triple Gaussian 


2" d order poly 


M 77 


double Gaussian 


Eq.© 


T(3S) - 


■> 7r + 7r~77^ + ^~ 


AM, 


triple Gaussian 


2 nd order poly 




Eq.© 


2 nd order poly 


T{3S) - 


-»■ 77M~V~ 


am; 


double Gaussian 


l s *order poly 


M 77 


double Gaussian 


Gaussian 



TABLE III: Efficiencies (e se i) and number of signal events 
(N) for each channel; upper limit at 90% CL (Nul) is given 
in parentheses. Uncertainties are statistical only. 



Transition 


Final state 


e se ; (%) 


A r 


r(2S) - 






39.1 


170061±413 


T(2S) - 


-> 77T(1S) 


7r + 7r _ 77^i + /i _ 
77M + M~ 


18.5 
37.2 


22±5 
90±14 


T(35) - 


-> 7r + 7r _ r(lS) 


7T + 7r _ e + e _ 


25.0 
42.8 


31330±186 
58500±247 


r(35) - 


-> ijr(is) 


7r + 7r _ 77e + e _ 

7T + 7r~77/i + /Lt - 

77M + M~ 


18.1 
8.9 
18.5 


4±2 (<8) 
4±2 (<8) 
7±11 (<26) 




A (GeV/cf ) 




7U7iyy 



(GeV/cf) 




> 20 



10 



A M; (GeV/c 2 ) 



0.5 0.55 0.6 

(GeV/c 2 ) 



(c) 




0.3 



0.35 



5 30 

2 20 
o 

^ 10 

60 

I 





(d) 







0.5 



0.6 



A M; (GeV/c z ) a 



M~ (GeV/c0 



am; 



and M 77 distributions for (a,b) Y(2S) 



FIG. 3: 

?7T(1S) — >■ 77/i + ^i _ candidates, and (c,d) T(3S) 
T]T(1S) — ► 77^i + /i~ candidates. Data are represented by dots, 
the fit results as solid curves and the background components 
by the dashed curves. 



> 

s 



1 



6^ 

4- 

2- 




0.25 0.3 0.35 0.4 0.45 

A H, (GeV/c 2 ) 




A (GeV/c 2 ) 



6 

4- 

2- 

L 



4- 



2- 




(d) 



0.5 0.55 0.6 0.65 

M^CGeV/c 2 ) 



(f) 



0.5 0.55 0.6 0.65 



M^CGeW) 



FIG. 2: AM, and M^ 77 distributions for (a,b) Y(2S) 
r/T(lS) — > 7r + 7r _ 77/i + pi~ candidates, (c,d) T(3S) 
V T(1S) 
V T(1S) 

by dots, the fit results as solid curves and the background 
components by the dashed curves. 



7r 7r 77/i (J, candidates, and (e,f) T(3S I ) — > 
7r + 7r _ 77e + e~ candidates. Data are represented 



tions. The systematic uncertainties on the AM WW and 
E* 12 vetoes for cross-feed dipion and radiative transi- 
tions are estimated by comparing the corresponding effi- 
ciencies in data and MC samples. In order to take into 
account possible discrepancies between simulation and 
data, the dipion events are generated using values for 
the transition matrix elements varied of ±ler with re- 
spect to those measured by CLEO [l(| ■ The difference in 
the efficiency is treated as a systematic uncertainty. The 
systematic uncertainty due to the choice of signal and 
background PDFs is estimated by using different func- 
tions, or different values for the fixed parameters. The 
complete list of contributions to the systematic uncer- 
tainty is summarized in Table IIV1 The total systematic 
uncertainty for each dataset is estimated by summing all 
the contributions in quadrature. 

The value of the branching fraction (B) , or upper limit 
on the branching fraction, for each mode is: 



B 



N 



(2) 



^sel X ^secondary 

where TV is the signal yield or upper limit on the sig- 
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TABLE IV: Sources of systematic uncertainty on the branching fractions B and on the ratios of partial widths, for each channel 
analyzed. All errors are given in percent. When both of the leptonic T(1S) decays are analyzed, the values in parentheses refer 
to the corresponding e + e~ final states. 





Y(2S) ->■ 


r(35) 




n + n~Y(lS) 






7r + 7r-T(lS') 






Source 






r l ~^ 77 






r/ — > 77 


Nr 




0.9 






1.0 




Tracking 


1.4 


1.4 


1.0 


2.5 


2.5 


1.7 






3.6 


3.6 




3.6 


3.6 


Lepton identification 




1.1 




1.0 (1.2) 


1.0 (1.2) 


1.0 


■k + tv~ model 


0.5 






0.4 (1.5) 






Selection 


0.4 


2.6 


5.5 


0.9 (1.2) 


4.4 (5.3) 


5.6 


PDFs 


0.1 


5.4 


5.0 


0.1 


5.4 


5.0 


Total B 


2.9 


7.6 


8.7 


3.6 (4.1) 


8.6 (9.1) 


8.1 


Total ratio 




7.2 


8.3 




8.3 (8.9) 


7.8 



TABLE V: Measured branching fractions and ratios of partial widths for hadronic T(nS) transitions. The first uncertainty is 
statistical, the second systematic. All ULs are at 90% of CL. The PDG values and the relevant predictions are given also. 





This work 


PDG [11] 


Predictions 5} 


B[Y(2S] -> TflTO-S)] (10" 4 ) 
B[Y(2S) -> 7r+7r-T(15)] (10~ 2 ) 
r[T(25) -> V Y(lS)} 3 

r[r(2S) -> n+n-r{is)] 1 ' 


2.39 ±0.31 ±0.14 
17.80 ±0.05 ±0.37 

1.35 ±0.17 ±0.08 


9 i +0.8 
z - L -0.7 

18.1 ±0.4 
1.2 ±0.4 


7-16 
40 

1.7-3.8 


B[Y(3S) -> ?7T(15)] (10" 4 ) 
B[T(3S) ^Tv+TY-r(lS)] (10~ 2 ) 
T[Y(3S) -> „T(LS)] 3 
r[T(35) -> n+n-T(lS)] 1 ' 


< 1.0 

4.32 ±0.07 ±0.13 

< 2.3 


< 1.8 
4.40 ±0.10 

< 4.2 


5-10 

5 

11-20 



nal yield. For a given channel, when both the leptonic 
T(IS) decays are available, their signal yields are first 
combined in a weighted average, where the weight is the 
inverse of the squared sum of the statistical and the sys- 
tematic uncertainties on each yield, considering only the 
systematic contributions that are uncorrelated. We as- 
sume B[T(1S) -> e+e - ] = B[T{IS) -> E3- For 
the rj transitions, the signal yields extracted from the two 
different rj decays are combined with the same weighted 
average technique. The results are shown in Table IVl 

We can also provide improved measurements of the 
differences between the Y invariant masses, using the fit- 
ted value of AM™ for both the Y(2S) -> tt+tt-TXIS') 
and the Y(3S) ir + ir~Y(lS) transitions. The val- 
ues are 562.170 ± 0.007(stat.) ± 0.088(syst.) MeV/c 2 and 
893.813 ±0.015(stat.)±0.107(syst.)MeV/c 2 , respectively, 
where the latter value is obtained as a weighted average 
of the values for the electron and muon samples. The 
systematic uncertainties are due mainly to the track mo- 
mentum measurement, which is related to the knowledge 
of the amount of detector material and of the magnetic 
field 0. 

We have presented a study of Y(3S) — > T(1S) and 
Y(2S) — > Y(1S) hadronic transitions. We have re- 



ported an improved measurement of B[Y(2S) — > r/Y(lS)] 
and a 90% CL UL on B[Y(3S) -)• r)T(lS)} compatible 
with, and more precise than, earlier measurements [2J, 
thus, further constraining theoretical predictions (see Ta- 
ble [Vj. We have also presented new measurements of 
B[Y{nS) -)• ir+TT-r(lS)} with n = 3,2, which we find 
to be compatible with earlier measurements [ill ]. Us- 
ing the independent BABAR measurement of B[Y(3S) — > 
XY(2S)] x B[T(2S) -)• 7r+7r-r(l,S')] in the inclusive di- 
pion spectrum [14{ . we extract the value B[Y(3S) — > 
XY(2S)} = (10.0 ± 0.6)%. 

Improved measurements of the ratios T[Y(nS) — > 
TjY(lS)]/T[r(nS) -)• 7r+7r-T(15)], for which systematic 
uncertainties partially cancel, have been presented also 
[llT |. The suppression of the Y(nS) — > rjY(lS) transi- 
tions with respect to the Y(nS) ir + ir~Y(lS) ones is 
confirmed to be higher than predicted by the QCDME 
[l[ and not compatible with other models 0, HI . 
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